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ABSTRACT 



Context. As the favoured progenitors of long-duration gamma-ray bursts, massive stars may represent our best signposts of individual 
objects in the early Universe, but special conditions seem required to make these bursters. These are thought to originate from the 
progenitor's rapid rotation and associated asymmetry. 
\ Aims. To obtain empirical constraints on the interplay between stellar rotation and wind asymmetry, we perform linear Hcf spec- 

• , tropolarimetry on a sample of 18 spectroscopically peculiar massive O stars, including OVz, Of?p, Oe, and Onfp stars, supplemented 

i-C ' by an earlier sample of 20 O supergiants of Harries et al., yielding a total number of 38 0-type stars. Our study's global aim is to 

characterize the differences between, and similarities amongst, different classes of peculiar O stars and to establish in how far they 
differ from garden- variety O stars. 

^ . Methods. Our linear (Stokes QU) spectropolarimetry data should be regarded a geometric counterpart to (Stokes /) spectral classifi- 

^ ' cation, setting the stage for circular (Stokes V) polarimetric searches for magnetic fields. 

^ ' Results. Despite their rapid rotation (with vsin/ up to ~400 km/s) most O-type stars are found to be spherically symmetric, but with 

notable exceptions amongst specific object classes. We divide the peculiar O stars into four distinct categories: Group I includes the 
suspected young zero-age main sequence OVz stars and related weak-winds objects, of which the magnetic star O' Ori C is the most 
famous member These objects show no evidence for significant linear polarization. Group II includes the Of?p stars, in which one 
^ ' of its members, HD 191612, was also found to be magnetic. These objects show more linear polarization activity than those in group 

00 I L Group III includes the Oe stars, which have been suggested to be the more massive counterparts to classical Be stars, and Group 

00 ■ IV concerns the Onfp stars. Objects from the latter two groups are on the high-end tail of the O star rotation distribution and have in 

OO ' the past been claimed to be embedded in disks. Here we report the detection of a classical depolarization "line effect" in the Oe star 

' HD453I4, but the overall incidence of line effects amongst Oe stars is significantly lower (1 out of 6) than amongst Be stars. The 

chance that the Oe and Be datasets are drawn from the same parent population is negligible (with 95% confidence). This implies there 
is as yet no evidence for a disk hypothesis in Oe stars, providing relevant constraints on the physical mechanism that is responsible 
for the Be phenomenon. Finally, we find that 3 out of 4 of the group IV Onfp stars show evidence for complex polarization effects, 
which are likely related to rapid rotation, and we speculate on the evolutionary links to B[e] stars. 
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X, 

H 1. Introduction scuration of the O star(-disk) system by their natal dusty cocoons 

- - I at their earliest times (e.g. Zinnecker & Yorke 2007). To estab- 

Massive stars have a pronounced effect on their surrounding in- jj^^ the evolutionary links in the early lives of O stars from their 

terstellar medium and on galaxy evolution, but despite their key prg.main sequence, to the zero-age main-sequence (ZAMS), to 

role in our Universe, huge uncertainties remain with respect to p^^^ main-sequence phases, one needs to characterize the differ- 

their formation and early evolution. As long as direct imaging ^^^.^^ between and similarities amongst their class, 
of the inner regions around young massive stars is beyond the 

capabilities of current imaging instrumentation, the tool of spec- O star spectral classification has been performed very suc- 
tropolarimetry may prove the most suitable avenue for obtaining cessfuUy over the last four decades. Important classes in this re- 
information about asymmetric geometries - of relevance if we spect are the so-called "Oe" and "Onfp" (Walborn 1973; also re- 
want to know which massive stars make long-duration gamma- ferred to as "Oef" Conti & Frost 1974) stars. The double-peaked 
ray bursts (GRBs). Balmer emission-line profiles of Oe stars and the He ii /14686 
For intermediate-mass young Herbig Ae/Be stars. Ha spec- emission-line profile of Onfp stars have been interpreted as sig- 
tropolarimetry has been very successful showing that Herbig Be "^"mg lom^ed disks, suggesting they might be higher mass ana- 
stars with masses up to -10-15 are embedded in disks on the log^es of classical Be stars (Conti & Leep 1974). Evidence for 
smallest spatial scales (Vink et al. 2002). For the more massive *e disk hypothesis has however, even after four decades, not yet 
O-type stars, the study of their youthful phase is more challeng- been found. 

ing due to the shorter evolutionary timescales as well as the ob- Almost two decades ago, Walborn & Parker (1992) iden- 

tified a group of potentially young O stars in the Magellanic 

Send Oj^nTj/ re^Me.sfj fo: lorick S. Vink, jsv@arm.ac.uk Clouds, the so-called OVz stars, which may be located close 
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to the ZAMS (hence the "z" designation in their spectral type). 
Walborn (2009) more recently presented a list of suspected 
ZAMS O stars, identifying three categories of ZAMS candi- 
dates: (a) a group of Ovz stars which show particularly strong 
He II /14686 absorption, indicative of a weak wind and low lumi- 
nosity (Walborn & Blades 1997), (b) stars with broad and strong 
hydrogen lines possibly indicating a high gravity, and (c) stars 
with very weak UV wind lines for their spectral types, indica- 
tive of a weaker wind due to a lower luminosity (see the lUE 
atlas of Walborn et al. 1985; Heydari-Malayeri et al. 2002). We 
note there is a difference between the phenomenology of "weak 
wind stars' which are objects with very weak UV wind signa- 
tures/or their spectral type and stars that are subject to the weak 
wind problem, which refers to an issue which seems to appear 
when comparing empirical and theoretical mass-loss rates for 
O-type stars below a luminosity of log (L/Lq) ^5.2 (Martins et 
al. 2005b, Mokiem et al. 2007, Puis et al. 2008, Marcolino et al. 
2009). 

Next to spectral Stokes I classification. Ho; spectropolarime- 
try adds geometrical information, as the linear polarized contri- 
bution of starlight has a more selective origin than total light. 
In its simplest form, the technique is based on the expectation 
that Ha line photons arise over a larger volume than stellar con- 
tinuum photons, such that line photons undergo fewer scatter- 
ings off the circumstellar disk than continuum photons, and the 
emission-line flux becomes less polarized than the continuum. 
This results in a smooth polarization variation across the Ha line 
profile: the "line effect". The high incidence of these line effects 
amongst classical Be stars (Vink et al. 2002 counted 26 out of 44 
from Poeckert & Marlborough 1976, i.e. ~55%) indicated that 
the envelopes of classical Be stars are not spherically symmetric. 
These findings are now taken as compelling evidence that clas- 
sical Be stars are embedded in circumstellar disks (e.g. Waters 
& Marlborough 1992, Quirrenbach et al. 1997, Oudmaijer et al. 
2005). 

In a similar vein, to characterize the early lives of massive 
O stars Ha spectropolarimetry may provide an appropriate tech- 
nique to explore source asymmetries, and to search the imme- 
diate ionized regions for disk signatures. Harries et al. (2002) 
performed Ha spectropolarimetry on a sample of 20 mostly nor- 
mal O supergiants, and found that in the majority of their sam- 
ple line effects were absent (75%), indicating that despite their 
rapid rotation the environments around O-type stars are gener- 
ally spherically symmetric. Here we present a linear Ha spec- 
tropolarimetry survey of 18 peculiar O-type stars. This supple- 
ments the sample of Harries et al. which consisted of 16 normal 
and 4 peculiar Onfp and Of?p stars, yielding a total sample size 
of 38, including 16 normal and 22 peculiar O-type stars. 

We divide the peculiar objects into four groups, (i) the sus- 
pected young zero-age main sequence (ZAMS) stars (OVz and 
possibly related objects), (ii) the peculiar Of?p stars (with opti- 
cal carbon iii emission lines as strong as the nitrogen iii lines), 
for which one of its members HD 191612 was recently found to 
be magnetic (Donati et al. 2006). (iii) the Oe stars - the suppos- 
edly hotter counterparts of classical Be stars, and finally: (iv) the 
rapidly rotating Onfp stars ("nfp" roughly indicating rapid rota- 
tion "n", strong wind emission "f" and spectral peculiarity "p")- 

The goal of our study is to establish the relationships be- 
tween the four groups of peculiar objects and to compare their 
linear polarization properties to those of normal O-type stars. In 
particular, we study the circumstellar environments for signs of 
youth and circumstellar geometries. Making these distinctions 
is expected to lead to a better understanding of both the early 
evolutionary path of massive stars towards the main sequence. 



as well as the subsequent rotational evolution of O stars beyond 
the ZAMS. This is particularly relevant with respect to the GRB 
phenomenon as special circumstances seem required to make 
them (e.g. Podsiadlowski et al. 2004). These are likely associated 
with the progenitor's rapid rotation, which could be obtained in 
both single (Yoon & Langer 2005, Hirschi et al. 2005, Woosley 
& Heger 2006) and binary scenarios (e.g. Cantiello et al. 2007, 
Wolf & Podsiadlowski 2007). The direct progenitors of GRBs 
are believed to be Wolf-Rayet stars, for which asymmetry sig- 
natures have been found in -15-20% of them using the tool of 
linear spectropolarimetry (Harries et al. 1998, Vink 2007). This 
might be accounted for with a scenario in which just the most 
rapidly rotating WR stars produce asymmetries. It seems only 
natural to assume that their rapidly rotating precursor O stars 
might show line effects too. 

The paper is organized as follows. In Sect. |2]we briefly de- 
scribe the observations, data reduction, and analysis of the linear 
polarization data, followed by a description of the resulting Ha 
line Stokes / and QU linear polarization profiles for the 4 dif- 
ferent subgroups of peculiar O stars (in Sect. O. In Sect. |4] we 
discuss the constraints these observations provide on the circum- 
stellar environments of the four groups of O stars. We summarize 
in Sect. |5] 

2. Observations, data reduction, and metliodology 

The linear spectropolarimetry data were obtained during two 
runs in Feb and Aug 2003 using the RGO spectrograph of 
the 4-metre Anglo-Australian Telescope (AAT). Two additional 
stars were observed using the ISIS spectrograph on the 4-metre 
William Herschel telescope (WHT) in 2003. The observations 
and data reduction of the AAT data is extensively described in 
Davies et al. (2005) whilst the same is done for the WHT data 
in Vink et al. (2005b). For the young star 0' Ori C the data are 
supplemented by three sets of WHT data from 1995/1996 (see 
Table 1), for which the data reduction procedure is described in 
Oudmaijer & Drew (1999). 

To analyze the linearly polarized component in the spectra, 
the spectrographs were equipped with the appropriate polariza- 
tion optics, consisting of a rotating half-wave plate and a cal- 
cite block to rotate and separate the light into two perpendicu- 
larly polarized light waves. For each exposure, four spectra are 
recorded: the ordinary (O) and extra-ordinary (E) rays of both 
the target and the sky. One complete observation set consists of 
a series of four exposures at half-wave plate position angles of 
0°, 45°, 22.5°, and 67.5° to obtain the linear Stokes parameters 
Q and U. Polarization and zero-polarization standards were ob- 
served regularly, revealing an intrinsic instrumental polarization 
of the order of 0.1 per cent. We did not attempt to correct for 
this, as our main aim is to investigate Hor polarimetric signatures 
relative to the continuum. 

The E and O ray data were imported into the polarimetry 
package ccd2pol, incorporated in the figaro software package 
(starlink). The Stokes parameters Q and U were determined, 
leading to the percentage linear polarization P and its position 
angle in the following way: 

p = V(e' + u^) (1) 

= i arctan(-|) (2) 

We note that a position angle (PA) of 0°, i.e. North-South on 
the sky, is represented by a vector that lies parallel to the posi- 
tive Q axis, whereas a PA of 90° (i.e. East- West) is positioned 
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in the negative Q direction. Positive and negative U axes thus 
correspond to position angles of respectively 45°and 135°. 

The data were subsequently analyzed using polmap. The 
achieved (relative) accuracy of the polarization data is in prin- 
ciple only limited by photon-statistics and can be very small 
(typically 0.01 %). However, the quality and the amount of data 
taken on spectropolarimetric standard stars is at present not yet 
sufficient to reach absolute accuracies below 0.1% (Tinbergen & 
Rutten 1997). 

We note that a non-detection would imply the wind is spher- 
ically symmetric on the sky (to within the detection limit). The 
detection limit is inversely dependent on the signal-to-noise ra- 
tio (SNR) of the spectrum, and the contrast of the emission line 
as the line-emission is depolarizing the (polarized) flux from the 
continuum. The detection limit for the maximum intrinsic polar- 
ization APiimit can be represented by: 



APii„,i,(%) 



100 



lie 



SNR lie - 1 



(3) 



where lie refers to the line-to-continuum contrast. This detec- 
tion limit is most useful for objects with strong emission lines, 
such as Ha emission in Luminous Blue Variables (Davies et al. 
2005), where the emission completely overwhelms underlying 
photospheric absorption. However, for the O-type stars studied 
in this paper, the line emission is much weaker (whilst some 
Group I objects even show \ia in absorption) and the detection 
limit equation, as provided by Eq.[3l looses its meaning. 

In general, we aim for an SNR in the continuum of 1000, 
or more, corresponding to changes in the amount of linear po- 
larization of 0.1%, or less. This way, we should be able to infer 
asymmetry degrees in the form of equator/pole density ratios, 
Peq/Ppoie of ~1.25, or larger (e.g. Harries et al. 1998), with some 
small additional dependence on the shape and inclination of the 
disk. 

Just as we do not correct for instrumental polarization, no 
correction for interstellar polarization (ISP) is made either, as 
the ISP only adds a wavelength-independent polarization vector 
to all observed line and continuum wavelengths. 

3. Results 

The linear continuum polarization of O-type stars is thought to 
be caused by scattering of stellar photons off electrons in the 
circumstellar environment, but in addition there may be an in- 
terstellar component to the measured level of polarization. The 
continuum (excluding the Ha line) polarization data for all our 
targets are summarized in Table [T] in the form of the mean R- 
band percentage polarization and its position angle 6 (columns 
4 and 5). Polarization variability is commonplace amongst O- 
type stars (e.g. Hayes 1975, Lupie & Nordsieck 1987), and ac- 
cordingly we do not look for perfect agreement between our 
continuum polarization measurements and those in earlier litera- 
ture. Nevertheless, the values of the polarization quantities %Pol 
and 9 are generally consistent with previous continuum measure- 
ments (see columns 6 and 7). 

Following Harries et al. (2002), we constructed maps of the 
ISP through searching the polarization catalogue of Mathewson 
et al. (1978) for objects lying within 5° from our targets. In 
some cases, the ISP maps appear both well-ordered and well- 
populated and in these cases column (10) states "OK". In case 
of sightlines where the ISP maps appear disordered (or sparsely 
populated) column (10) states "No". In those cases where the 
match between the measured polarization and the ISP is good. 



the overall polarization is assumed to be dominated by the ISP, 
and the intrinsic polarization will be small. 

The Ha line shapes of our 0-star targets include absorption, 
emission, double-peaked emission, and P Cygni profiles. Plots of 
the spectropolarimetric data are presented in the different panels 
of Figs.[T]to|6] The polarization spectra are presented as triplots, 
consisting of Stokes I, P, and 9. 

3.1. OVz stars 

The membership criteria for OVz and related ZAMS candidates 
have been summarized in Sect.[T]and Walborn (2009). 

The Ho- line polarimetry data for the 9 ZAMS candidate O- 
type stars are very similar in the sense that Ha is mostly in ab- 
sorption (see Fig|2]i. In none of them do we find evidence for a 
variation in the polarization percentage or position angle across 
the line. In order to address the question of whether this implies 
the objects are intrinsically unpolarized, we compare the contin- 
uum polarization data to previous measurements, as well as to 
the ISP 

We start off with the two Orion ZAMS candidate objects 
©' Ori C and 0^ Ori A. As these two objects are members of 
the Orion nebula cluster (ONC) star forming region, it should 
be no surprise that the ISP maps for the objects appear disor- 
dered. The same is true for the four OVz stars in the Carina star 
forming region, c.f. HD 93128, HD93129B, CPD-58°261 1, and 
CPD-58°2620, but for these four objects we have no information 
regarding previous linear polarization data. 

In our AAT data of 0' Ori C and Ori A, the Ha lines 
are in absorption, and although this does not necessarily imply 
there cannot be any line effects, it is still the case that line effects 
are easier to detect across Ha when the line is seen in emis- 
sion, which happens on a regular basis for the oblique rotator 0' 
Ori C (Stahl et al. 1996, Donati et al. 2002, Gagne et al. 2005). 
Multiple observations of 0' Ori C over four different epochs are 
plotted in Fig.[T] The Stokes I Ha line is seen to vary from one 
night to the next from a pure emission line on 31 Dec 1995 to an 
inverse P Cygni profile on 1 Jan 1996. In our data from 31 Dec 
1996 and 14 Feb 2003 the line displays a pure absorption pro- 
file. We note that the very narrow features that appear close to 
line centre are simply the result of imperfect nebular subtraction. 
Despite the large differences in the character of the Ha line, the 
linear polarization properties of 0' Ori C (see Table 1) do not 
appear to vary significantly. 

Given that our continuum polarization measurement for the 
other Orion object Ori A appears to have changed slightly 
with respect to earlier measurements (column 7), this might in- 
dicate that the spectroscopic binary ©- Ori A is variable in its 
intrinsic linear polarization properties, which might be consis- 
tent with wind-wind interaction or polarization variations due to 
a magnetically confined wind, such as in the Bp object crOri E 
(Townsend & Cohen, in prep.). 

For the remaining three candidate ZAMS O-type stars, cf. 
HD 42088, HD 54662, and HD 152590, our PA values are con- 
sistent with earlier measurements and the ISP, and it is probably 
safe to state that the measured polarization for these objects is 
predominately of interstellar origin, with only an exceedingly 
small intrinsic polarization component. 

3.2. Of?p Stars 

The mysterious Of?p phenomenon is referred to when the carbon 
m optical emission at ~4650 A is equally strong as that of the 
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Table 1. Polarization data for four groups of peculiar O-type stars observed during the AAT and WHT runs described in the body 
text. 



Name 


Spec.Tp 


Date 






P'" (%) 


gilt 


P'^ (%) 


7^ 


ISP map 


(I) OVz and weak wind dwarfs: 




















0' Ori C 


07 Vp 


31-12-95 


0.31 ± 0.01 


139.9 ± 0.3 
















01-01-96 


0.29 ± 0.01 


148.2 ± 0.8 
















31-12-96 


0.29 + 0.01 


144.0 ± 0.6 
















14-02-03 


0.284 ± 0.005 


136.3 ± 0.5 


0.35 


151.8 


0.240 


88.5 


No 


0^ Ori A 


09.5 Vpe 


14-02-03 


1.051 +0.005 


99.7 ±0.1 


0.702 


93.1 


0.221 


88.5 


No 


HD 42088 


06.5 V 


14-02-03 


2.431 ± 0.008 


177.0 ± 0.1 


2.53 


178 


2.126 


168.8 


OK 


HD 54662 


06.5 V 


14-02-03 


0.881 ± 0.098 


145.5 ± 3.8 


0.83 


145 


0.481 


154.4 


No 


HD 93128 


03.5 V((f+)) 


12-02-03 


2.265 + 0.014 


89.0 ± 0.2 












HD93129B 


03.5 V((f+)) 


12-02-03 


1.889 + 0.063 


117.5 ± 1.0 












HD 152590 


07.5 V 


14-02-03 


0.998 ± 0.009 


29.1 ±0.3 


0.96 


19.2 


0.737 


27.9 


No 


CPD-58°2611 


06 V((f)) 


13-02-03 


2.222 ± 0.010 


84.7 ±0.1 












CPD-58°2620 


06.5 V((f)) 


14-02-03 


1.717 ± 0.009 


90.6 ±0.1 













(11) Of?p: 



HD 148937 



06.5 f?p 



12-02-03 
15-08-03 



1.580: 

1.803 : 



0.020 
0.014 



46.9: 

49.3 : 



0.4 

:0.2 



1.63 



50.6 



1.727 



36.1 



OK 



(III) Oe stars: 



HD 45314 


09:pe 


14-02-03 


1.316± 


0.091 


167.4 


± 2.0 


1.38 


171 


1.369 


171.6 


OK 


HD 46485 


07 Vn(e) 


WHT 03 


1.887 ± 


0.013 


26.1 


±0.2 


1.89 


159.0 


0.808 


165.8 


OK 


HD 60848 


08 V:pev 


WHT 03 


0.431 ± 


0.007 


135.7 


± 0.4 


0.78 


56 


0.372 


27.0 


No 


HD 120678 


08 Illmep 


13-02-03 


1.637 ± 


0.015 


77.8 


±0.3 


1.7 


76 


1.781 


71.7 


OK 


Zeta Oph 


09.5 Vnn 


14-02-03 


1.572 ± 


0.004 


127.3 


±0.1 


1.42 


127 


0.411 


76.4 


No 


HD 155806 


07.5 V[n]e 


14-08-03 


0.821 ± 


0.013 


156.4 


± 0.5 


0.86 


147 


0.742 


159.0 


OK 



(IV) Onfp: 



HD 152248 



07Ib:(n)(f)p 15-08-03 0.477 ± 0.016 114.9 ±1.0 



0.66 



112 0.812 



27.9 



OK 



The spectral types are taken from the O star catalogue of Maiz-Apellaniz et al. (2004) and listed in column 2, with the dates of the observations 
provided in the third column. The errors in the polarization data (column 4) are of the order of 0.01% based on photon-statistics only. Yet, 
systematic (external) errors in the polarization are estimated to be 0.1%. The errors in the Position Angle 6 (column 5) are of order a degree. The 
literature values of the polarization and PA are taken from the catalogue of Mathewson et al. (1978) and so are the values for the ISP (columns 
8 and 9). The last column (10) indicates whether an ISP map is well ordered and well populated (see body text in Sect. |3]for more detailed 
explanation). 



nitrogen iii emission at 4640A (see e.g. Maiz-Apellaniz et al. 
2004). Only three such objects are known in the Milky Way. 
These are HD 191612, which is magnetic (Donati et al. 2006), 
HD 108, and HD 148937. 

Linear spectropolarimetry data on two different epochs in 
Feb and Aug 2003 for HD 148937 are presented in Fig.|3] along- 
side data for HD 108 from Harries et al. (2002). Although both 
epochs of HD 148937 show a P Cygni profile in Ha, there are 
notable changes in the absorptive part of the P Cygni profile. 
Neither epoch however shows a significant change in the polar- 
ization character across the Ho-. The continuum PA of the two 
datasets has varied only slightly, from = 49.9° + 0.4° in Feb 
2003 to 6» = 49.3° + 0.2° in Aug 2003, i.e. the same within the er- 
rors, although the level of polarization may have increased from 
P = 1.58% to 1.80% over the same time interval, but as the sys- 



tematic eiTor for these kind of data is of order 0. 1 % (see Sect.|2]l, 
this only constitutes a 2cr effect. 

3.3. Oe stars 

Oe stars are defined as near main sequence O-type stars with 
their Balmer lines in emission, possibly representing the higher 
mass analogues of classical Be stars (Walborn 1973, Negueruela 
et al. 2004). 

The spectropolarimetry data for the six Oe stars are depicted 
in Fig. m The Stokes I normal intensity data show a variety 
of line profiles, involving absorption, single, and double-peaked 
emission lines. HD45314 shows a clear polarization variation 
across the Stokes I emission profile that is as broad as that of 
the emission line itself. We can safely consider this effect to be 
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1 


^ \ — 


— ' — \ — 










Fig. 1. Polarization spectra of the Group i peculiar O star 0^ Ori C at four differentepochs. Stokes I spectra are shown in the lowest 
panels of the triplots, the levels of %Pol in the middle panel, whilst the PAs {0; see Eq. 2) are plotted in the upper panels. The 
data are rebinned such that the Icr error in the polarization corresponds to 0.05% or 0. 1 % as calculated from photon statistics. The 
narrow features that appear in three of the Stokes / line profiles are due to imperfect nebular subtraction. 



hd 45314 



o 




0.8 



1 .2 



q (%) 



Fig. 5.QU diagram of the Group m Oe star HD 453 14. Note that 
the normalized Stokes parameters u - U/I and q - Q/I are 
depicted. The continuum polarization appears on the top of the 
diagram at u ~-0.45%. The Ha line excursion is shown (with an 
arrow) downwards to a value of u as low as — 0.9%. 



the result of the classical "depolarization" line effect (see the 
extensive discussion in Vink et al. 2002). In such a case, one 
would expect to observe a change in %Pol across the line that 
follows the shape of the intensity profile, but with the PA re- 
maining constant. However as a result of the vector addition of 
interstellar polarization, it may be observable in the PA instead. 
Indeed, although the line effect is present in both the polarization 
percentage (middle panel) and the PA (upper panel), the varia- 
tion in the PA is not intrinsic to the star but a result of inter- 
vening ISP, which simply shifts the datapoints across the QU 



plane. This conclusion is confirmed and explained when we al- 
ternatively plot the same data in a QU diagram (see Fig|5]). In 
this representation the continuum polarization appears at the top 
of the diagram, corresponding to u ^-0.45%, whilst the Ha line 
makes a linear excursion downwards (with the arrow denoting 
the wavelength direction) to a value of u of ^-0.9%. In other 
words, the intrinisic polarization of HD 453 14 is of order 0.45% 
at an intrinsic PA of ^135°. 

The other five Oe stars are not subject to any polarization line 
effects across Ha, although their Stokes I profiles show a variety 
of shapes. This is an intriguing result, esp. when one considers 
the fact that most vsin/ values are large (see Table 2) and double 
peaked emission is seen in non-line-effect objects like HD 60848 
and HD 155806, which have strong Stokes I Ha emission lines. 
In the past, double-peaked emission lines in Stokes I intensity 
have often been interpreted as signalling circumstellar disks, but 
here we are reminded by the fact that double-peaked profiles 
can for example also be produced in an expanding shell that is 
spherically symmetric and would, irrespective of the inclination 
of the stellar rotation axis, not produce any linear polarization. 

With respect to the continuum polarization measurements 
listed in Table 1, we count three objects for which we trust the 
ISP estimates. For the case of HD 46485 the measured PA ap- 
pears different from prior measurements, as well as from the ISP. 
For HD 120678 and HD 155806 there is no notable difference 
between our PA measurements and earlier literature, nor do we 
note any variation with respect to their ISP estimates. We can 
safely assume these two objects are unpolarized. For the case of 
HD 60848 our PA measurement of 136° deviates notably from 
a previous measurement of 56° but as the level of polarization 
is low, the error in the PA is quite large. The continuum PA 
measurement of 127° for ^Oph is consistent with the previous 
measurement, however the continuum polarization percentage 
has dropped to 0.6%, compared to 1.4% on a previous occasion, 
leaving open the possibility that this object is intrinsically polar- 
ized but that the wind emission at the epoch of observation was 
too weak to produce a line effect. 

We summarize the Oe star polarization results in the follow- 
ing way. Out of the sample of six Oe stars, we can state that at 
least two of them are upolarized. If one were willing to accept 
time variability or ISP corrected polarization results, the remain- 
ing four Oe stars from our sample may be intrinsically polarized. 
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Fig. 2. Polarization spectra of the remaining Group i ZAMS candidate O stars. Note that the spectrum of HD93129B is dominated 
by its 1.5 mag brighter companion HD93129A at only 3 arcsec away. 



hd 148937 12/02/2003 



hd 148937 15/08/2003 




6560 6580 
Wavelength (1) 



Fig. 3. Polarization spectra of the Group n Of?p stars HD 148937 (during both the Feb and Aug 2003 runs) and HD 108 (from 
Harries et al. 2002). The data are rebirmed such that the Icr error in the polarization corresponds to 0.05% or 0.1% as calculated 
from photon statistics. 
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Fig. 4. Polarization spectra of the Group m Oe stars. 



However, when we do restrict ourselves to the original criteria 
for detecting intrinsically polarized O-type stars, only accepting 
objects that show polarization changes across spectral lines as 
sufficient evidence for intrinsic polarization, we find that only 
one out of six has a line effect. This fraction is statistically lower 
than that in classical Be stars, which have an incidence of line 
polarization effects of ~55% (Poeckert & Marlborough 1976, 
Vink et al. 2002). The implications of this result are discussed in 
Sect, a 

3.4. Onfp stars 

The Onfp stars are defined as objects in which He n 4686A 
shows centrally reversed emission (Walborn 1973), sometimes 
referred to as Oef stars (Conti & Frost 1974). In general these 
objects rotate rapidly, although there might be a variety of phys- 
ical phenomena that give rise to the strong He ii emission, such 
as a strong stellar wind, or gas streams in close binaries. 

The spectropolarimetry data for the Onfp star HD 152248 is 
presented in Fig.|6] alongside data for three additional Onfp stars 
from Harries et al. (2002). HD 152248 shows the kind of P Cygni 
Stokes I Ha profile that are reminiscent of the more famous Onfp 
supergiants A Cep and ( Pup. The difference between them is that 
the very high signal-to-noise data for A Cep and ^ Pup show line 
effects (see Sect. |4] for a discussion on the nature of these line 
effects), whilst the lower signal-to-noise data of our HD 152248 
data only shows a tentative indication of a small variation in the 
PA, but the signal is not high enough to count HD 152248 as a 
line-effect star at the present time. 



4. Discussion 

In the following we place our results in the broader context of 
relevant characteristic information known about the four sub- 
groups of O-type stars. This might provide constraints on for- 
mation models for massive stars, as well as the relevance of stel- 
lar rotation and magnetic fields with respect to disk formation 
around massive stars. 

To help assessing the probability of finding a line effect in 
an individual object or a subgroup of peculiar O-type stars, we 
quantify a representative wind density of our objects, pwind = 
M/{R^Voo) (column 9 of Table |2|l. Note that we choose theoreti- 
cal over empirical mass-loss rates (as determined from quantita- 
tive spectroscopy), as only 5 out of these 23 targets have empir- 
ical mass-loss rates available derived from line-blanketed mod- 
els (see Repolust et al. 2004, Martins et al. 2005b, Mokiem et 
al. 2005, Marcolino et al. 2009). Furthermore, the rates derived 
from different spectral ranges appear to be highly uncertain. For 
instance, for ^ Oph, Marcolino et al. (2009) provide a mass- 
loss rate as low as log(M) = -8.80 from the ultraviolet, whilst 
Mokiem et al. (2005) give log(M) = -6.84, determined from Ha. 
The latter value is in very good agreement with the theoretical 
value of log(M) = -6.83 (see column 8). In other words, the use 
of theoretical rates is not only preferred for practical reasons. We 
note that if the terminal velocity is not available, we multiplied 
the escape velocity by the canonical factor 2.6 for stars on the 
hot side of the bi-stability jump (Lamers et al. 1995). 

As can be noted from column 9 of Table |2l the average 
wind density per subgroup does not appear to be very differ- 
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ent. Interestingly, the line-effect Oe star HD45314 is actually 
an object with one of the smallest wind densities given in col- 
umn 9. We can thus conclude that the average wind density does 
not play a major role in determining whether or not a certain 
subgroup of O stars is subject to linear polarization. Rather, it 
is predominately determined by the equator/pole density ratio, 
Peq/Ppoie which is required to be larger than 1.25, in order to 
achieve the measurable amount of linear polarization of 0. 1 % or 
more (see Sect. 2). 

4. 1 . Group I: OVz and weak-wind stars 

Given that this subgroup may contain a number of young O- 
type stars, one might possibly expect that some could be sur- 
rounded by circumstellar material leftover from the star for- 
mation process - in a manner analogues to pre-main sequence 
Herbig Ae/Be stars. For this reason, we searched the 2MASS 
catalogue for JHK magnitudes, but we did not find OVz stars to 
have a large near-infrared excesses (of more than ~ 0.5 mags), 
from which we conclude that it is unlikely that a large fraction of 
OVz and related objects be surrounded by large quantities of cir- 
cumstellar dust. However, gaseous circumstellar disks could po- 
tentially be detected using linear spectropolarimetry, if (i) these 
objects are sufficiently young, and (ii) if massive O stars form 
via disk accretion. 

A second reason why one might anticipate the presence of 
a gaseous circumstellar disk in these objects is that its most 
famous member 0' Ori C has a strong dipolar magnetic field 
(Donati et al. 2002, Wade et al. 2006) and probably a magnet- 
ically confined wind leading to a pile-up of gaseous material 
around the magnetic equator like in the Bp star cr Ori E (e.g. 
Babel & Montmerle 1997, ud-Doula & Owocki 2002). 

We can summarize the results of our linear Ha line spec- 
tropolarimetry of the nine Group i OVz and related stars with 
the simple statement that none of them shows strong evidence 
for linear polarization, and we consider a number of different 
reasons why this might be so. 

Perhaps the simplest explanation would be that the absorp- 
tive character of the Stokes I Ha profile prevents our technique 
from being sensitive enough. However, it would not explain the 
non-detection in the Dec 1995 data of ©' Ori C when Ho' was 
seen in emission (see Fig.[Tll- Even if our linear polarimetry tool 



fails to detect disks when Ha is mostly photospheric, and even if 
accretion disks are absent, it is still somewhat puzzling that even 
for 0' Ori C in which the presence of a magnetically confined 
wind seems well-established, we do not detect any evidence for 
flattened circumstellar material in our OVz hnear polarimetry 
data. 

It is possible that the circumstellar disk structures are dis- 
rupted by the magnetic pressure of the central object and that 
gaseous material is only present at larger radii, with an inner 
hole present closer to the star, where the number of free elec- 
trons would be lower than for normal O stars, limiting the level 
of linear continuum polarization. In other words, the absence of 
line effects in ©' Ori C might actually be consistent with a mag- 
netically confined wind scenario, where the magnetic pressure 
is too high close to the surface, creating an inner disk hole. It 
would be interesting to provide further constraints on the size of 
such a hole and this might yet again be achieved using line po- 
larimetry. Differences in the shapes of intrinsic line polarimetry 
(in contrast to line depolarization where the continuum is polar- 
ized) involving either single or double loops in the QU plane 
have the potential to constrain the size of disk inner holes (Vink 
et al. 2005a). This could work in case the line under considera- 
tion is in emission, which occurs quite regularly for objects like 
©' Ori C (see Fig. [TJ as well as the Bp star cr Ori E. It will be 
interesting to see whether this technique could be applied to the 
bulk of the Group i OVz and related objects. 

We note that Smith & Fullerton (2005) re-examined the ve- 
locity behaviour of UV wind lines in 0' Ori C and found that 
in contrast to cr Ori E, the UV profiles of 0' Ori C show sur- 
prisingly little evidence for the presence of a pile-up of material, 
suggesting that despite the strong dipolar magnetic field of ©' 
Ori C, it is too early to conclude the star is simply a massive 
analog of Bp stars. To continue the discussion along this line 
of reasoning, the small variations in the continuum polarization 
of ©2 Ori A may be related to the complex gas flows around 
the object. We recommend the Group i OVz and related objects 
be monitored simultaneously in linear and circular polarization 
to obtain a more comprehensive picture of the geometry of the 
magnetic field in conjunction with that of the circumstellar ma- 
terial. 
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Name log(L/Lo) T^if logg M R v^o logM p„ind vsin; 
(kK) (Me) (R^) (km/s) (M^yr') (IQ-" Mpyr-'/km s'' V) (tan/s) 



Group I OVz and weak wind dwarfs: 



0' Ori C 


5.13 


39 


4.1 


30 


8 


580 


-5.7 


530 


24 


0- Ori A 


4.62 


30.5 


3.92 


17 


7 


2500' 


-7.59 


1.9 




HD 42088 


5.23 


38 


4.0 


32 


10 


1900 


-6.19 


37 


60 


HD 54662 


5.2 


36.8 


3.92 


29 


10 


2456 


-6.40 


17 


73 


HD 152590 


5.00 


34.4 


3.92 


24 


9 


1785 


-6.64 


23 


60 


HD 93128 


5.66 


46.5 


4.00 


40 


10 


3210' 


-5.53 


100 


100 


HD93129B 


5.86 


44.0 


3.92 


52 


13 


3210' 


-5.25 


110 




CPD-58°2611 


5.3 


38.2 


3.92 


32 


10 


2870' 


-6.23 


21 




CPD-58°2620 


5.2 


36.8 


3.92 


29 


24 


1760' 


-6.43 


15 


110 



Group II Of?p: 



HD 148937 


5.75 


41.0 


4.0 


-55 


15 


2285 


-5.45 


71 


76 


HD 108 


5.40 


37.0 


3.75 


-35 


12 


1960 


-5.97 


37 


78 


HD 191612 


5.4 


35 


3.5 


-35 


15 


2450' 


-6.14 


17 





Group III Oe stars: 



HD 39680 


5.5 


40 




45 


11 


1675 


-5.69 


90 


218 


HD 45314 


4.72 


31.5 


3.92 


18 


8 


2410' 


-7.37 


3.0 


285 


HD 46485 


5.1 


35.5 


3.92 


27 


9 


1800 


-6.44 


25 




HD 60848 


4.9 


33.4 


3.92 


22 


8.5 


1765 


-6.84 


12 


240 


HD 120678 


4.9 


33.4 


3.92 


22 


8.5 


2580' 


-7.02 


5.5 




Zeta Oph 


4.88 


32.1 


3.62 


20 


9 


1505 


-6.83 


12 


372 


HD 155806 


4.9 


33.4 


3.92 


22 


8.5 


2460 


-7.02 


5.5 


91 



Group IV Onfp: 



HD 152248A 


5.69 


33.3 


3.4 


41 


21.1 


2420 


-6.10 


16 


159 


HD 210839 


5.83 


36.0 


3.55 


62 


21.1 


2250 


-5.47 


33 


200 


Zeta Pup 


5.90 


39.0 


3.60 


54 


19.4 


2250 


-5.24 


68 


219 


Cyg OB2 No.5 


5.69 


33.3 


3.4 


41 


21.1 


2240' 


-5.64 


27 





The stellar parameters are taken from the following papers: Repolust et al. (2004), Martins et al. (2005b), Puis et al. (2006). For all stars without 
recent line-blanketed models, the calibrations of Martins et al. (2005a) are used. The mass-loss rates are from Monte Carlo radiation-driven wind 
models (Vink et al. 2000). The wind density is given as p„jnd = M/(R^v^). Where spectral subtypes are not provided, we interpolate, vsin; are 
mostly taken from Howarth et al. (1997). 

' Note that the terminal velocity is calculated by multiplication of the escape velocity by the canonical factor of 2.6 (Lamers et al. 1995). 



4.2. Of?p Stars 

The Of?p stars appear to be relatively slow rotating stars (see 
Table 2) that show dramatic periodic spectral variability (e.g. 
Naze et al. 2008). Contrary to the OVz object 0- Ori A, the 
PA of the Of?p star HD 148937 is constant with time, although 
it may be subject to variations in the level of linear polarization. 
Out of the 20 supergiants studied by Harries et al. (2002), the 
one Of?p star that was studied, HD 108, was one of only five 
objects to show a line effect. The character of the HD 108 line 
effect is consistent with classical depolarization (see Fig.O 

The third known Galactic Of?p star, HD 191612, has not yet 
been observed in linear spectropolarimetry mode, but interest- 



ingly it is one out of only three O-type star£] in which a mag- 
netic field has been discovered through circular spectropolarime- 
try (Donati et al. 2006, Howarth et al. 2007). Similarly to the 
comments made above with respect to the Group i objects like 
0' Ori C, it would be most informative to monitor the Of?p stars 
simultaneously in linear and circular polarized light. The fact 
that there is no evidence for rapid rotation in any of the known 
Of?p stars (see Table |2]i might be consistent with a scenario of 
magnetic braking for these objects. 



' The other two O-type stars with confirmed magnetic fields are 0' 
Ori C (Donati et al. 2002, Wade et al. 2006), and f Ori (Bouret et al. 
2008b) 
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4.3. Oe Stars 

The Oe stars have been suggested to be related to classical Be 
stars. The definition of classical Be stars is that they are near- 
main sequence B stars (of luminosity class m - v) in which the 
hydrogen Balmer lines are or have been in emission. The emis- 
sion lines often show blue- and red-shifted peaks, consistent with 
the presence of a rotating circumstellar disk (Struve 1931), and 
this disk model was confirmed decades later through both polar- 
ization (e.g. Poeckert & Marlborough 1976, Wood et al. 1997) 
and interferometry measurements (Quirrenbach et al. 1997). The 
situation for the Oe stars is however less clear On the basis of the 
Stokes 1 spectrum, Conti & Leep (1974) suggested these objects 
simply represent the higher mass analogues of the classical Be 
stars, but in more recent work of Negueruela et al. (2004) it was 
found that many of these objects were in fact spectroscopically 
classified too early, with the possibility that the more massive 
counterparts to classical Be stars, i.e. classical Oe stars, simply 
do not exist. 

The relevance of the general absence of Oe stars is that it 
could provide insight into the physical mechanism underlying 
the Be phenomenon, as some theoretical models for explaining 
the Be phenomenon predict disks to be present across the en- 
tire spectral range, whilst other models predict the presence of 
disks solely at certain spectral types, corresponding to specific 
effective temperatures, and possibly restricted to specific stellar 
masses and/or ages (e.g. Fabregat & Torejon 2000, McSwain & 
Gies 2005, Martayan et al. 2006). Reasons for an absence of the 
hotter analogues of classical Be stars could involve the more in- 
tense ionizing radiation or stronger radiatively-driven wind of 
the hotter O stars, weaker wind braking for B stars (e.g. Maeder 
& Meynet 2000), or special circumstances, such as wind bi- 
stability (Pauldrach & Puis 1990, Vink et al. 1999) to explain 
the presence of circumstellar disks exclusively amongst objects 
that fall in the B-star spectral range. 

Even though Negueruela et al. (2004) found some of the Oe 
stars to be classified "too early", it is still true that some of them 
fall within in the O-star spectral range, and as a population they 
are rapid rotators (see the large vsin/ values in Table|2]i. The rele- 
vant question is now whether these objects are embedded in geo- 
metrically thin disks, analogous to those present around classical 
Be stars. When we take our result of one out of six line depolar- 
ization effects in the Group iii Oe stars at face value and compare 
this incidence to that amongst the classical Be stars, whilst prop- 
erly accounting for the different sample sizes, we find that the 
chance the line effect data of the Be and Oe samples are drawn 
from the same parent sample is negligible, at 95% confidence. 
In other words, contrary to common belief, there is no evidence 
that early-type Oe stars such as HD 155806 are embedded in 
circumstellar disks. 

We subscribe to the view of Negueruela et al. (2004) that the 
group of Oe stars, and especially those with spectral types earlier 
than 09.5, deserve more detailed investigations to provide better 
constraints on the physical origin of the Be phenomenon. 



4.4. Onfp stars 

With respect to the incidence of line polarization effects for the 
Onfp stars, the situation appears to be different from that in the 
Oe stars discussed above. Out of the four Onfp stars studied 



spectropolarimetricalljQ, at least three of them show evidence 
for line effects across Ha. 

We should note that the physical phenomena that give rise 
to the peculiar He ii emission might be due to differing physical 
origins. Cyg OB2 No. 5 is classified as 07 lanfp, but this ob- 
ject does not display the characteristic Onfp He ii 4686A profile. 
Instead, it is a peculiar and variable interacting binary. The He ii 
4686 A emission in HD 152248 may also be due to gas streams 
in a close binary. This object was recently studied by Mayer et 
al. (2008) and the stellar mass of ^30 Mq in that paper is prob- 
ably more accurate than the stellar mass listed in Table 2 (as it 
is simply derived from the binary motion). However, we find the 
Hff shape in HD 152248 very reminiscent of that in A Cep and ( 
Pup. 

With respect to the nature of the line effects in Onfp stars 
we cannot draw any definitive conclusions, yet we make the fol- 
lowing remarks. The fact that the polarization changes in A Cep 
occur predominantly across blue-shifted absorption might indi- 
cate that it is not the line depolarization effect that is at work 
here. Instead, these polarization profiles seem to be more remi- 
niscent of the "McLean" effect (McLean 1979) in Herbig Ae/Be 
stars (Vink et al. 2002) or of optically pumped gas (Kuhn et al. 
2007). However, the most convincing modelling so far suggests 
that these line effects are the result of the star's rapid rotation, 
which causes an asymmetry in velocity space (Harries 2000). 

Given the small sample of Onfp stars studied so far and 
the complexities involved in the line polarization profiles of the 
Group IV Onfp stars, it would be premature to conclude that 
the current sample of line polarization effects need necessar- 
ily indicate the presence of disks around Onfp stars. It should 
also be noted that geometrical effects other than disks, such as 
wind clumping (see Davies et al. 2007 and references therein) 
might also give rise to the observed levels of linear polariza- 
tion. Indeed, it may also present an alternative explanation for 
the peculiar Stokes / spectrum of ^ Pup (Bouret et al. 2008a). 
Nonetheless, similarly to Group iii Oe stars. Group iv Onfp stars 
appear to be rapid rotators (see Table |2]i, with the difference be- 
ing that the incidence of line effects in Onfp stars appears to 
be consistent with that expected for objects surrounded by disks 
with random orientations, although this is based on extremely 
low number statistics. 

Given their rapid rotation and the range of their intrinsic lu- 
minosities (log(L/L0) ^5.7 - 5.8; see Table 2), it is tempting 
to speculate the Galactic Onfp stars are the precursors of the 
Galactic B[e] supergiants, which are also believed to be rotat- 
ing rapidly but for which the evidence of equatorial disks seems 
better established (Zickgraf et al. 1985). Interestingly, the B[e] 
phenomenon has been found to extend to significantly lower lu- 
minosities in the Large Magellanic Cloud (Gummersbach et al. 
1995) and it will be informative to see if the Onfp phenomenon 
is also observed towards lower luminosities in the LMC (see 
Walborn 2008). 

Finally, it will be crucial to enlarge the database of linear 
spectropolarimetry data for these Group iv Onfp stars, as they 
appear to be the most promising subgroup of O-type stars of be- 
ing a significant disk producing massive star population, with 
potential implications for constraining progenitor models for 
gamma-ray bursts and other exotic phenomena that may be re- 
lated to rapid rotation. 



- Note that we have now counted ^ Pup as an Onfp star. According to 
Walborn's spectra the object is not Onfp, but Conti noticed the charac- 
teristic He II 4686A emission in his Oef notation, indicating variability. 
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5. Summary 

We have presented linear spectropolarimetry data on a sample 
of 18 peculiar O-type stars, supplementing an earlier sample 
of 20 mostly normal O stars of Harries et al. (2002), yielding 
a total of 38 O-type stars. On the basis of the spectral Stokes 
I peculiarities we divided these peculiar objects into four 
separate groups and studied their linear Stokes QU polarization 
properties. 

Group I included the suspected young zero-age main sequence 
OVz stars and related objects. For this group we can summarize 
the results with the simple statement that none of them shows 
evidence for significant amounts of linear polarization. It is not 
inconceivable that circumstellar disks would be disrupted by the 
magnetic pressure of the central object, and the absence of polar- 
ization line effects might be considered consistent with a mag- 
netic scenario, where the magnetic pressure could be responsible 
for an inner hole, which limits the number of free electrons close 
to the stellar surface and the associated level of linear continuum 
polarization. 

We also suggested that the measured variations in the 
continuum polarization of objects like 0^ Ori A may be related 
to the complex gas flows associated with its binary nature or the 
formation of magnetically confined disks and we recommend 
Group I objects be simultaneously monitored in linear and 
circular polarized light. 

Group II included the spectroscopically variable Of?p stars, 
which seem to comprise a group of slowly spinning stars that 
are subject to dramatic periodic spectral variability. 

Contrary to the Group i stars like 0^ Ori A, the PA of the 
Of?p object HD 148937 is constant with time, although it may 
be subject to variations in the level of linear polarization (a 2 cr 
effect). Interestingly, a second member of the group HD 108 was 
one of the few O-type stars in Harries et al. that showed a line 
effect. The third member HD 191612 was not observed in linear 
spectropolarimetry mode, but interestingly it was recently found 
to be magnetic, possibly providing a link to some of the younger 
magnetic stars from group i such as ©' Ori C. 

Given the high incidence of both linear and circular polar- 
ization plus the fact that there is no evidence for rapid rotation in 
any of the Of?p stars (see Table |2]i, leads us to speculate that it 
is the presence of a magnetic field that is the underlying reason 
for this peculiar spectral classification. 

Group III includes Oe stars, which have been suggested to be 
the more massive counterparts to classical Be stars. However, 
when we take our linear Ho- polarimetry results at face value 
and compare the low incidence (one out of six) with that of the 
classical Be stars (26 out of 44), thereby properly accounting for 
the different sample sizes, we conclude that the chance that the 
line effect Oe and Be samples are drawn from the same parent 
sample is negligible at the 95% confidence level. 

In other words, there is as yet no evidence that stars with 
spectral types different from the B-type range possess circum- 
stellar disks. 

Group IV concern the rapidly rotating Onfp stars. Similarly to 
objects from group III, evidence for a disk hypothesis has yet 
to be found for this subgroup. By contrast to the Oe stars, out 
of four Onfp stars studied, at least three show evidence of line 
polarization effects. Due to small number statistics it is not pos- 
sible to conclude that our results imply the presence of disks. 



but interestingly the high incidence of line effect Onfp stars does 
appear to be in line with that expected for objects with random 
orientations (Poeckert & Marlborough 1976, Vink et al. 2002). 

We speculated that the Galactic Onfp supergiants could be 
the precursors of B[e] supergiants, given their similar luminosity 
range (around logiL/Lo) ^5.7-5.8), their rapid rotation, and 
possibly the presence of circumstellar outflowing disks in both 
object classes. It will be interesting to enlarge the current linear 
polarimetry dataset for group iv Onfp objects, as these objects 
appear to be the most promising subgroup of O-type stars with 
line effects. 

To conclude, we presented the first linear QU polarization 
data for O-type stars that exhibit Stokes / peculiarity. We showed 
that those objects which are peculiar in Stokes / have a higher 
chance of showing QU polarization. Future Stokes V circular 
polarimetry studies of both normal and peculiar O-type stars e.g. 
via the mimes collaboration with CFHT/Espadons should teach 
us how the incidence of magnetic fields compares to the Stokes 
/ and Stokes QU phenomenology described here, as the rota- 
tional evolution of massive O stars is probably driven by a com- 
plex interplay between mass loss, rotation, and magnetic fields. 
This interplay needs to be constrained empirically and under- 
stood theoretically before we are able to explain which massive 
stars produce exotic phenomena such as long-duration gamma- 
ray bursts. 
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